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Recently, the strain state of a piezoelectric electrode has been found to impact the electrochemical activity 
taking place between the piezoelectric material and its solution environment. This effect, dubbed 
piezocatalysis, is prominent in piezoelectric materials because the strain state and electronic state of these 
materials are strongly coupled. Herein we develop a general theoretical analysis of the piezocatalysis process 
utilizing well-established piezoelectric, semiconductor, molecular orbital and electrochemistry frameworks. 
The analysis shows good agreement with experimental results, reproducing the time-dependent voltage 
drop and H 2 production behaviors of an oscillating piezoelectric Pb(Mg 1/3 Nb 2 /3)03-32PbTi03 (PMN-PT) 
cantilever in deionized water environment. This study provides general guidance for future experiments 
utilizing different piezoelectric materials, such as ZnO, BaTi0 3 , PbTi0 3 , and PMN-PT. Our analysis 
indicates a high piezoelectric coupling coefficient and a low electrical conductivity are desired for enabling 
high electrochemical activity; whereas electrical permittivity must be optimized to balance piezoelectric and 
capacitive effects. 

Piezocatalysis is a new approach toward enabling or enhancing electrochemical processes by making use of 
the strain state of a piezoelectric material 1 . Piezocatalysis is the product of an intimate interaction between 
the native electronic state of the piezoelectric material, the chemistry of the surrounding medium, and a 
strain induced piezoelectric potential. The action of mechanically deforming a piezoelectric material induces a 
perfuse electric field which augments the energetics of both free and bound charges throughout the material 1,2 . 
The thermodynamic feasibility and kinetics of electrochemical processes occurring at the surface of the piezo- 
electric material sensitively depend upon the electrochemical potential difference between charges on the piezo- 
electric's surface and in the surrounding medium 3 6 . Thus piezoelectric potential, which can dramatically affect 
the difference between these electrochemical potentials, is a new means of modulating the material's electro- 
chemical activity via its strain state. 

Recently, a piezocatalysis process was demonstrated by our study of a strained ferroelectric Pb(Mg 1/3 Nb 2 / 3 )0 3 - 
32PbTi0 3 (PMN-PT) beam in a deionized (DI) water system, from which a strong dependence of hydrogen 
evolution from the aqueous surroundings on the material's piezoelectric potential was observed 1 . In addition, 
numerous recent works have confirmed the correlation between electrochemical activity and piezoelectric or 
ferroelectric polarization in a broader sense. For example, a study conducted using ferroelectric poly(vinylidene 
fluoride) (PVDF) has demonstrate that in-situ piezopotential can influence lithium battery charging behavior 7 . 
Electrochemical deposition was found to be selectively activated by the ferroelectric domain polarization 8-15 . 
However, to date there is no general theoretical analysis of the piezopotential's effect on electrochemical activities, 
e.g. how one piezoelectric material's activity differs from another; the influence of metallic electrodes as compared 
to bare piezoelectric surface on reaction output, and how free charge in piezoelectric material systems affects the 
piezocatalysis process. In this paper, we address the piezocatalysis system in generality to elucidate the details 
underlying these open questions and illuminate trends for further experimental study. 

In order to clearly understand the piezocatalysis process, a conventional electrocatalysis process is discussed 
first, where the application of electrical potential from an external power source is a typical means of driving the 
electron transfer reactions (Fig. la). The applied potential can result in one of the following two processes: (1) 
lowering electronic energy levels of unoccupied states within the electrode to a magnitude less than that of the 
highest occupied molecular orbital (HOMO) in solution; (2) raising occupied states within the electrode above the 
lowest unoccupied molecular orbital (LUMO) in solution. Under the first condition, electrons will leave 
the HOMOs in solution and transfer to the unoccupied states within the electrode - oxidizing the solution 
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Figure 1 | Energy diagrams describing the electrochemistry and piezocatalysis process, (a) The effect of applying a sufficient positive or negative bias to 
an electrode is to increase anodic and cathodic current, respectively, (b) In the case of a bare piezoelectric, conduction and valence bands act as the 
reservoirs for electrons donated or accepted from molecules in solution. The piezoelectric polarization (P v -£) applies a variable bias across the material, 
lifting and lowering valence band and conduction band energies, (c) Applying electrodes to the piezoelectric simplifies the electron reservoir to that of the 
metal's Fermi energy, while maintaining the piezoelectric potential as the source of bias. 



(center panel of Fig. la). Under the second condition, electrons will 
leave the occupied states in the electrode and transfer to the LUMOs 
in solution - reducing the solution (right panel of Fig. la). When the 
electrode is a non-metallic solid, the e(p HOMO and e<l> LUM0 logic 
extends to the materials valence (e(9 VB ) and conduction band 
(e?) CB ) edges, respectively. 

For piezocatalysis, the external power source is replaced by piezo- 
electric potential which results from the piezoelectric polarization 
(Ppz). A perfectly insulating piezoelectric material is the most ideal 
and simplest case and it will be analyzed first. For such a material, 
mechanical deformation creates a perfuse electric field inducing a 
total energy shift across the material, V p , given by: 



V„ 



w x T k d 



<-xk 



(1) 



where is an applied stress in the k dimension, d x/ t are the piezo- 
electric moduli, So is the electrical permittivity of free space, and is 
the relative permittivity in the x dimension and W x is the width of the 
piezoelectric material in the x dimension. In a one-dimensional case, 
x and k are equal and the subscripts are dropped. 

This piezoelectric potential, which changes the energetics of the 
valence band (VB) and conduction band (CB) across the piezo- 
electric material, can have a dramatic effect on the material's inter- 
action with its environment (Fig. lb). In this model it is assumed that 
stain does not change the magnitude of the band gap. If the e<p vs 
approaches e(j) LUMO , it becomes energetically favorable for electrons 
to leave the VB and enter the LUMO (center panel of Fig. lb). If 
the etp cs approaches e#> HOMO , it becomes energetically favorable for 



electrons to leave the HOMO and enter the CB (right panel of 
Fig. lb). Placing metal electrodes between the piezoelectric and solu- 
tion, with their continuous density of states about their Femi ener- 
gies, simplifies this situation: the piezopotential now acts as a bias, 
lifting and lowering the metal's Fermi energy e<p M (Fig. lc). 

Piezoelectric materials, which are the source of bias in the case of 
piezocatalysis, are capable of achieving extremely high potentials 
(tens to hundreds of volts) when subjected to moderate to severe 
strain 16 . Under such circumstances (i.e. where electrode potentials 
versus standard hydrogen electrode (SHE) exceeding ~3 volts) many 
chemical species in contact with the piezoelectric will be thermody- 
namically capable of undergoing reduction or oxidation reactions. 

In order to deduce the maximum quantity of oxidation or reduc- 
tion reactions possible for a given deformation, it is necessary to 
quantify the interaction between charge exchange at the piezo- 
electric/electrolyte interface and the effect of charge exchange on 
the piezopotential. A characteristic of piezocatalytic systems is the 
piezopotential drop that takes place during charge transfer to and 
from its surfaces (or electrodes). In this way the system acts as a 
capacitor. The rate of piezopotential drop is dependent upon both 
the piezoelectric material's properties and the nature of the solution 
in which it is submerged. In the case where the concentration of 
reactive species in solution is low, e.g. when pure Milli-Q water is 
used for the hydrogen reduction reaction, and the power available for 
driving electron-transfer reactions is sufficiently high, the rate of 
electrochemical reactions at the electrode is dominated by the dif- 
fusion rate of redox species to the electrode. In the case of Milli-Q 
water, the redox species are protons, hydroxide and impurity species. 
Mass transfer to the electrode is described by the Nerst-Planck 
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equation, which in its one-dimensional form along the x-axis is 
written as: 

0Q(x) 

C'X 



Ji(x)=-Di- 



— DiCi(x) + CiV(x) 



where Jj(x) is the flux of species i at distance x from the surface, D; is 
SCi(x) 

the diffusion coefficient, — - — is the concentration gradient at 

ox 

distance x, F is the Faraday constant, _R is the gas constant, T is 

d(p(x) 

absolute temperature, — - — is the potential gradient, Q is the con- 
centration of species i, and v(x) is the velocity with which a volume 
element in solution moves along the axis 17 . The expression describes 
the contributions of diffusion, migration and convention, respect- 
ively, to the flux of species i. 

Neglecting migration and convective phenomena, the flux of spe- 
cies i is dependent upon diffusion alone. In the diffusion limited 
regime, the maximum rate of electron transfer from the piezoelectric 
is equal to the rate of reactant diffusion to the piezoelectric, resulting 
in a current density j given by: 

ji(x)= -niFDj— - — (3) 



where nj is the number of electrons per reaction event with species i". 

Applying appropriate boundary conditions to equation (3) and 
taking account of the capacitive nature of the piezoelectric bias 
source, yields an expression of the piezopotential V p as a function 
of time (a variant on the Cottrell equation): 



SVp 
8t 



v.F^niD^qfi 



(4) 



where C; is the bulk concentration of species is a parameter taking 
on a value between 0 and 1 that toggles the kinetics of the electrodes' 
reactivates with species i, and t is time 18 . 

Under conditions of low reactant and electrolyte concentrations 
and high positive electrode potential (e0 i(/Alo >2eV vs SHE), a kin- 
etics parameter (/j) less than 1 causes a semi-diffusion controlled 
regime to form where charged reactant species (e.g. protons and 
metal ions) capactively couple to the electrode's surface, effectively 
reducing the surface potential before they are electrochemically 
reduced or oxidized (Fig. 2a). Under the application of a large pos- 
itive electrode potential, a dilute (e.g. Milli-Q water) system's capa- 
citance is well approximated by the Helmholtz model 17 . The voltage 
drop expected in time t by both electron transfer reactions and capa- 
citive effects is given by integrating equation (4): 



2Ff'/ 2 

£(>Er,xTC 1/2 



v x nPt^cf, + w H Y, ztmD^cft. -fi)] (5) 

;'=/ i=l 

where Zj is the charge sign (1 for cation, — 1 for anion) of species i, w H 
is the thickness of the Helmholtz layer, and I is the number of dif- 
ferent species in solution. The positive potential of the electrode 
dictates that a positive charge in solution will contribute to a positive 
capacitive current, lowering the piezopotential. The first term in 
equation (5) describes the piezopotential change as redox reac- 
tions proceed, the second term describes the potential change due 
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Figure 2 | The voltage drop across a strained piezoelectric is dependent upon the rate of capacitive and Faradic effects, (a) The charged piezoelectric's 
surface induces both capacitive double layer effects and electron transfer events across the interface, (b) At high negative piezoelectric potentials, the 
electrode is thermodynamically capable of reducing all reported impurity ions in solution. As the potential increases below -2.0 volts, the majority of 
side reactions are no longer possible, (c) Depending on t p and the values of kinetic parameters f ;il and f i)2 , the rate of voltage drop experienced as a 
function of time can take on a wide range of values, (d) Applying a recursion procedure to the data obtained by Wang et al, values of t p = 0.042 s, 
fi,i = 0.715 and f; >2 = 0.07 are determined. 



SCIENTIFIC REPORTS | 3 : 2160 | DOI: 1 0. 1 038/srep02 1 60 



3 



to capacitive coupling of not yet reacted species at the piezoelectric 
material's surface. 

Equation (5) governs how the potential on the piezoelectric's wall 
should drop from the time (t = 0) of initial mechanical deformation, 
when the piezopotential can be a value of tens or hundreds of volts, to 
a time (f = t p ) when the potential has been reduced to a reasonable 
value about the SHE (e.g. ±2 Vvs SHE). Between time 0 < t< f p ,fhe 
electrode potential is sufficiently high as to reduce or oxidize species 
within solution in an effectively nonselective manner. At t > f p , the 
piezoelectric potential falls within the range of typical ei^ HOMO and 
&j>iUMo va l ues of ions in solution resulting in a more selective reduc- 
tion and oxidation process (Fig. 2b). 

The processes described above pertain to both surfaces of the 
piezoelectric material; on one surface oxidation dominates while 
reduction dominates on the other. The full voltage change, as mea- 
sured across the piezoelectric material, is a combination of both of 
these processes happening simultaneously. An electrode potential 
above (in the case of oxidation current) or below (reduction current) 
the selective potential window of ~ ± 2 V versus SHE here are treated 
identically. This results in a total voltage reduction across the piezo- 
electric material of twice the expected magnitude from a single wall 
(equation (5)). 

As the piezopotential reaches a modest range about 0 V (t = f p ), 
the electrochemical activity of some ions is effectively switched off 
while others will remain highly active. Those ions which are no 
longer electrochemically active are still free to participate in capacit- 
ive effects (Fig. 2b and second term in equation (5)). The voltage 



decrease across the piezoelectric from r„ 



oo is thus due to a 



separate set of electrochemical and capacitive effects conducted by a 
subset (i...m) of the total chemical species present (i. . . /). By applying 
equation (5) to both the selective and unselective regime we can 
describe the magnitude of piezopotential at all times by the following 
derived expression: 



V p = 



-H(t p -t) 



4R 1 / 2 

1/2 



Wh 



H(t-t p ) 
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i= i 



/= 1 
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£ «A 1/2 qfi,i + — £ z."A l/2 Ci(l ~fi,i) 
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i = l 



Sh 



Sh 



where H(x) is the Heaviside function defined as 0 when x is less than 
0, and 1 when x is equal to or greater than 0. is the electrical 
permittivity of the Helmholtz layer. The third and fourth terms in 
equation (6) acts to add continuity between the selective and unse- 
lective regimes. 

Equation (6) is a simplification of the activity occurring around 
the piezoelectric as the piezopotential drops to within the vicinity 
of e^ HOMO and e$ LUM0 of species in solution. As the energetic 
advantage for electron transfer diminishes, the current density j 
begins to show its exponential dependence upon driving force 



(e<p HOMO — s(j> LUM0 ). The current density exponentially depends 
upon applied potential about E e : 



j=jo exp 



-exp 



where j 0 is the exchange current density and is equal to: 

(l-a fl )nFE e 



./o : 



(7) 



(8) 



f p and the 



n is the number of electrons transferred in the reaction, k° a and k a c are 
the rate constants for the anodic and cathodic current, respectively, 
Cr and Co are the concentrations of the reduced and oxidized spe- 
cies, respectively, a a is the anodic transfer coefficient, E e is the redox 
potential of a species in solution, and if is the applied bias {i.e. Fermi 
energy difference between eq> m and £ c ) ls . For the purposes of simpli- 
fication, the effects of equation (7) will be approximated by utilizing 
two different kinetics constants acting in equation (6) for 1 > / > m: 
/u and f 2 ,i for the nonselective time interval 0 £ t 
selective interval f p £fat», respectively. 

To elucidate the behavior of the system, Fig. 2c is a plot of equation 
(6) as it is applied to the case of a piezoelectric material (PMN-PT) 
strained in Milli-Q water. We take a value of electrical permittivity 
Spmn-pt,x = 8000 and a thickness of PMN-PT w th = .23 mm 19 . The 
value of the Helmhotlz layer thickness (wjj) was fixed to 2.75A with a 
Eh = 80. Under standard conditions, Milli-Q water has a resistivity 
value of 18 fl - cm due to the contribution of / "impurities", including 
hydroxide and hydronium (c 12 = —0.100 ppb'mol'cm -3 , D 12 = 
9.3* 10~ 5 cm 2 sec), organics (c 3 = 1.106 ppb - mol - cm~ 3 , A =0.52* 
10 -5 cm 2 sec), silicates (c 4 = 0.5534 ppb - mol - cm~ 3 , A = 1.28* 
10 _5 cm 2 sec), and heavy metal ions (calcium, sodium, chloride, etc; 
c 5 = 0.1106 ppb-mol-cnr 3 , A = 1.2 * 10~ 5 cm 2 sec) 20 . Based on 
experimental measurements discussed in the following section, tp 
is chosen at 0.042 s and/ 1( is 0.715 (between 0>f>f p , 71.5% of 
all species immediately undergo redox reactions upon their contact 
with the electrode, 28.5% of species first capacitively couple to the 
electrode). 

The physical consequences of the voltage decrease described by 
equation (6) are: (1) the creation of reaction byproducts in the med- 
ium surrounding the piezoelectric material (e.g. hydrogen gas, oxy- 
gen gas, chlorine gas, oxidized organics, etc); (2) the buildup of a 
capacitive layer in solution alongside the piezoelectric material; and 
(3) a reduction in the electric field present inside the piezoelectric 
material. These effects will continue to change with time until the 
potential decrease (equation (6)) approaches the total potential ori- 
ginally generated (equation (1)). The rate of these effects depend 
sensitively upon the rate constants/^ and/ 2 ;, the concentration of 
reactive species and their diffusion coefficients. The natural driving 
force for these chemical changes is the piezopotential which, as 
described by equation (1), depends upon the electrical permittivity, 
mechanical strain, Young's moduli and piezoelectric moduli. This is 
true regardless of whether the electrochemically active medium is in 
contact with metallic electrodes or the piezoelectric material's bare 
surfaces. In the case of bare piezoelectric surfaces, the band gap and 
positions of e^ CB and e<t> VB relative to e^ HOMO and &j> LUM0 also 
determine the nature of what happens during straining. 

When comparing the piezocatalytic abilities of various piezo- 
electric materials, it is useful to choose a performance metric from 
equation (6) such as the quantity of a reaction byproduct, e.g. H 2 gas. 
This value depends sensitively upon the parameters we choose. 
Previous work reported the production of hydrogen gas via piezo- 
catalysis between a PMN-PT single crystal slab (2 mm X 10 mm in 
size) and Milli-Q water 1 . Using recursion we are able to fit equation 
(6) to the piezopotential verses time profile obtained from a working 
electrode on a PMN-PT surface under one strain (Fig. 2d). The 
results yield t p = 0.042 s,/ H = 0.715, and f 2i = 0.07 under the 
constraints that after t = t p only H + reduction and H 2 0 oxidation 
were possible at the electrodes, whilst unreacted H + , OH~, silicate 
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and heavy metal ions only acted capacitively. These low kinetics 
values likely resulted from a low concentration of electrolyte. The 
switch over at f p = 0.042 s from one reaction regime to the other 
occurred at a potential (difference) value between the electrodes of 
approximately 2 V. A gold electrode in equilibrium with pH 7 Milli- 
Q water is in equilibrium with the H + reduction and H 2 0 oxidation 
reactions. This places the Au electrode's potential (¥>ElectrodeHOMo) at 
0.41 Vvs. SHE (0.615 V below the ep LUMO offf+ and -0.615 above 
the ei^ HOMO of H 2 0). A symmetric voltage shift applied to both Au 
electrodes of 1 V at t p means an overpotential (^Electrode homo — 
^lumo) °f — 0.385 V for the H + reduction reaction and an over- 
potential (^Electrode homo ~ Plumo) of °- 3 85 V for the H 2 0 oxidation 
reaction. These are the potential values taken at which the system 
switches from a nonselective regime to a selective regime. 

In order to derive a function of H 2 evolution dependency on strain, 
the minimum strain S m i n ,i capable of driving electrochemical reac- 
tions for an arbitrary piezoelectric is necessary and determined by 
augmenting equation (1): 

c _ £ 0 £ r,x(yElectrodeHOMO ~ ^LUMo) , Q \ 

>->min,l — „, (y) 

2w x Yd xk 

where Y is the Young's modulus of the material. For metallic Au 
electrodes, the preceding analysis determined the potential necessary 
to drive electrochemical H 2 production under the observed condi- 
tions to be 0.615 (^EkctrodeHOMO-^LUMo)- In the case of a bare 
piezoelectric, ( Electrode homo -Plumo) becomes (^vb-^lumo)- 
S m in,2 is denoted as the strain necessary to transition from the select- 
ive to nonselective regime, e.g. where KkctrodeHOMO-^LUMo) in 
equation (9) becomes ( Electrode homo -?>lumo +0.385V). Using 
the value parameters fit from experiment in conjunction with the 
values of piezoelectric materials' parameters listed in Table 1, the H 2 
generation capacity per straining event (HMetaljotal) for a multitude 
of Au-electrode coated piezoelectric materials is depicted in Figure 3a 
and is given by a manipulation of equation (6): 

HMetal,Total=H(S — S m i n ^)HMetal,SmallsH(S m \ n ^—S) + 

(10) 

XargeS 

HMetal,Smatts is the H 2 production for piezopotentials that fall within 
the selective regime and is given by: 



H-Metali 



11 



where t x is the time required to exhaust the hydrogen production 
capability for a piezoelectric operating only within the selective 
regime (S min ,i <S<S min , 2 ) and is given by: 



ti = 



reiectrodeHOMO 



-fLUMO^OV^ 



n H + w x d]jI C h + fi,H+ Wh 



\ 



— + — £A 1/2 c,(W 2 ,,)+ £ A 1/2 c, 



(12) 



where e^o is the electrical permittivity of water. HMetaljMrgeS is the H 2 
production for piezopotentials that reach the unbiased regime 
(S>Smin,2) and is given by: 



Huetai 



LargeS ' 



(13) 



where f 2 (equation (14)) is the time it takes for the piezopotential to 
fall from an arbitrarily large potential (V pmax ) to that which brings it 
into the selective regime (i/'lumo"'" Vop' where (p 0f is the overpoten- 
tial necessary to transition between reaction regimes, e.g. 0.385 V), 
and f 3 (equation (15)) is the time it takes for that same piezoelectric to 
reduce its potential from 9>lumo + Vop t° 9>lumo- 



h = 



and 



<(<P: 



ElectrodeHOMO 



4F 



1 't,.v 



i—l ^H 2 o i—i 



<p Op £ 0 ^Tl 

IF 



(14) 



V 



E H 20 



^D," 2 Q(l-/ 2 , i )+ £ D] /2 Q 



J 

(15) 



At large strains (strain > .00025), PZT dominates H 2 gas production 
per unit strain with a value of 9.1 X 10 14 H 2 - strain, PMN-PT and 
BTO have similar production rates of 5.4 X 10 14 H 2 - strain and 4.5 X 
10 14 H 2 - strain respectively, while ZnO maintains a modest H 2 pro- 
duction rate of 79 X 10° H 2 - strain. From equation (13) it can be 
determined that the H 2 production depends upon the relative weight 
of/i,H+ to/2,H+> dxk and e r , x . When/ijj+ is much larger than/ 2 jj+, 
tJT-l dominates the expression where a large value results from a 
larger V p combined with a large e r , x . However e rjc is contained 
within the denominator of the expression determining V p (equa- 
tion (1)). Thus, while a large d^ is always valuable to the piezo- 
catalysis process, a balance must be achieved with e r x . This is why 
PZT, possessing both large and moderately large e r>x , has the 
largest H 2 production rate. Physically this competition represents a 
balancing between the piezoelectric, voltage generating capabilities 
of the material, and the capacitive nature of the material. The inset in 
Figure 3a enlarges the small strain region (strain < 0.00004), dem- 
onstrating the dynamics of initiation H 2 production. Materials with 

the smallest ratio, the most prominent being ZnO, begin H 2 

djd 

production at strains far smaller than a material like PMN-PT. 
This is because the thermodynamic accessibility of electrochemical 
reactions depends upon the relative energetics of donor (ef/? HO Mo) 
and acceptor (e(j> LUMO ) states and a high ratio of piezoelectric moduli 
to electrical permittivity strongly couples the strain state of the 
material to these state energetics (equation (9)). 

The predicted H 2 production capacity is compared to experimental 
data obtained from an oscillating gold coated PMN-PT cantilever 
with a peak piezopotential of 20 volts (Fig. 3b) 1 . At 20 Hz, each 
straining action takes place during 0.025 seconds, which is insuf- 
ficient for depleting all piezopotential and reaching a thermodynamic 
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Figure 3 | The H 2 production per straining event predicted as a function 

of strain magnitude, (a) In the case of piezoelectric materials with metal 

electrodes on their surfaces. Inset shows the turn-on strain value is 

dependent on the ratio of -y- . (b) Comparison between experimental data 
dj* 

of hydrogen generation over the course of an oscillating PMN-PT 
cantilever experiment and the predicted results from simulation under the 
conditions of poor kinetics (/[ = 0.07). Inset shows the envelope of 
hydrogen generation predicted from simulation under the case of good 
kinetics (fi = 0.715) and poor kinetics, (c) In the case of bare surfaced 
piezoelectric materials. Inset shows the turn-on strain values depends upon 
both the ratio of and the value of <p 0 „. 

equilibrium. Nonetheless, an approximation of the experiment can be 
constructed by truncating the nonselective regime at 0.025 seconds 
and applying a window of acceptable values on the kinetics parame- 
ter from 0.07 to 0.715 (inset of Fig. 3b). Experimental data fall 
around the lower bound of predicted H 2 evolution rate, exhibiting 
good agreement. This indicates that there remains much room for 



experimentally improving the H 2 production rate by means of adding 
a proton reduction co-catalyst and thus increasing the/j value. 

Several factors change when comparing a piezoelectric material 
coated with metal electrodes to a naked, insulating one. In general, a 
naked and perfectly insulating piezoelectric material does not neces- 
sarily possess a continuum of states about the redox potentials in 
solution. To achieve piezocatalysis (including H 2 production) with 
such a piezoelectric material, the conduction band and valence band 
must be moved sufficiently in potential so that their energies coincide 
with redox potentials under examination. Taking the production of 
H 2 as our prototypical case, the energetics of top-most valence band 
electrons (e<p VB ) must be lifted above the (e(f> LUM o) of H 2 . The min- 
imum strain, and thus potential, required to achieve this criteria is 
again given by equation (9), except in the case of a bare piezoelectric 
^ElectiodeHOMO becomes <p VR . The values of potentials necessary to 
cause the reduction of H + and oxidation of H 2 0 for various piezo- 
electric materials can be found in Table 1. The hydrogen production 
capacity of a bare piezoelectric material can be calculated by modi- 
fying equation (10): 



H Insulator, Total 



min,l 



- S)Hi mu l ator s m allsH(S — S m in,2) + 



where Hi„ su i ator Smalls is the H 2 production for piezopotentials that fall 
within the selective regime (S milU <S<S min ,2). H !nsu i atorMrseS is the 
H 2 production for piezopotentials capable to reaching the nonselec- 
tive regime (S>S mill , 2 ). 

Applying the values from Table 1 and the same kinetic parameters 
as in the metal electrode case (f u = 0.715, / 2i = 0.07, and (p 0p = 
— 0.385 V), Fig. 3c shows the H 2 production rate for various naked 
piezoelectric materials. Subjected to large strains (strain>.001), both 
naked and Au coated piezoelectric materials perform consonantly. 
Discrepancies between these two are more prominent under small 
strain (strain<0.0002), which are a result of the difference in the 
band structures of individual piezoelectric materials. Under small 
strain, each piezoelectric's H 2 -production's turn-on strain depends 
not only upon the — ^- ratio but also upon the value of the individual 

piezoelectric material's (eq> vs — etp LUMO ) value. 

Concerns over the stability of a piezocatalysis system comprising a 
bare piezoelectric material can be addressed by comparing piezo- 
catalysis with a controlled corrosion process, where reduction reac- 
tions take place through the piezoelectric's oxidation and vice versa. 
The maximum piezoelectric charge density at the surface is on the 
order of 10 13 charge'cm 2 , two orders of magnitude lower than the 
average atomic density of a solid-state surface. With oxidation events 
of this magnitude, it is expected that a piezoelectric capable of with- 
standing photocatalytic reactions with water is also able to survive 
the piezocatalysis process. 

In contrast to the previous analysis, most piezoelectric materials 
are not perfect insulators, mobile electrical charges inside the piezo- 
electric material respond to and rearrange according to its internal 
electric field. In bulk systems, these mobile charges act as extremely 
effective screening agents of the piezopotential, resulting in an effec- 
tive piezopotential given by a manipulation of the Gouy-Chapman 
capacitance model: 

2fcTsinh~ '[t^a&Y \/ 'e 0 e, r , x kn 0 T] 
2\/2s 0 £ rx kn 0 T 



Vse, 



where k is Boltzmann's constant,, and n 0 is the bulk concentration of 
free charge 17 . According to equation (17), a highly insulating bulk 
piezoelectric is required to reach an appreciable potential (Fig. 4a). 
Taking the values from Table 1 for various piezoelectric materials 
subjected to a strain of 0.002, the maximum concentration of free 
charges allowed in order for the piezoelectric materials' surfaces to 
obtain the H 2 production potential were 5.22 X 10" 40 , 2.33 X 10~ 27 , 
9.37 X 10~ 23 , and 1.38 X 10~ 26 m~ 3 , for ZNO, PMNPT, BTO and 
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Figure 4 | Free charge can dramatically affect the piezocatalytic ability of 
a strained piezoelectric. For a strain= .0.002, the effect of free charge on (a) 
piezoelectric potential and (b) H 2 production for various piezocatalysis 
systems. 



PZT respectively. The H 2 production capacity for these piezoelectric 
materials subjected to 0.2% strain was determined as a function of 
free charge by manipulating equation (16): 



Hsemi,Total — H (n maxl — n 0 )H Hi gi,„ 0 H(n 0 — « max ,2) + 
H(n 0 



(18) 



where Hmghn a is the expression for H 2 production while the piezo- 
potential is within the selective regime and Hi ow „ 0 describes the H 2 
production for a piezoelectric that has achieved a potential for acces- 
sing the nonselective regime. « max ,ii s the maximum number of free 
charges allowed for a given piezoelectric material under 0.2% strain 
to reach the threshold potential (^ V b — Plumo) f° r driving H 2 
production. n max ,2 is the charge concentration necessary to reach 
the potential threshold between the two reaction regimes 
[(^vb - ^lumo)+ Vop\- in equation (18) the dependence of H 2 pro- 
duction on n 0 comes exclusively from the fact that in the case of a 
semiconducting piezoelectric material, Vp max [equation (13) & (14)] 
becomes Vsemi [equation (17)]. 

The effect of free charge on H 2 production for various piezoelectric 
materials is demonstrated in Fig. 4b, where the hydrogen production 
for a strain of 0.2% is shown as a function of mobile charge concen- 
tration, n 0 . Behavior of particular note is the dramatic increase in H 2 
production that results from decreasing the mobile charge concen- 
tration ~5 orders of magnitude below the minimum concentration 
required to reach the H 2 potential for the piezoelectric materials. 
Additional mobile charge reduction beyond that initial 5 orders 
reduction has relatively little effect on H 2 evolution. 



Discussion 

A number of factors have been shown to augment the efficacy of 
strain to induce electrochemical reactions. The relative energies of 
states within the electrode with respect to HOMO and LUMO ener- 
gies in solution can both dramatically change the rate of chemical 
reactions and be a determining factor in which chemical reactions are 
allowed to proceed. These energy state positions depend directly on 
the magnitude of strain and piezoelectric coefficient while depending 
inversely upon the electrical permittivity. In the presence of free- 
charge, generated either by doping, photo or thermal excitations, 
the piezopotential can be markedly decreased with direct repercus- 
sions on reactivity. In order to increase the electrochemical activity of 
a strained piezoelectric, the value electrical permittivity needs to be 
optimized. The piezoelectric, semiconductor and molecular orbital 
frameworks discussed herein have historically and successfully been 
applied to solid state and electrochemical systems. It is expected that 
these frameworks can be extended to the piezocatalysis theory. The 
theoretical work presented here can provide fundamental guidance 
for designing and understanding the novel piezocatalysis phenom- 
enon from broad electrochemistry and piezoelectric material 
systems. 

Methods 

All calculations were performed in Wolfram Mathematica 8. To calculate the total H 2 
output from an insulating piezoelectric with metal electrodes, we first calculated the 
H 2 generated as the piezoelectric potential dropped from the maximum potential 
induced by strain [equation(l)] to the potential present at time t v . Then, the H 2 
generated after time t p was calculated using another value of the kinetic parameter. 
Summation of these two H 2 quantities gave the total H 2 output. The first time segment 
utilized kinetics parameter/i = 0.715; while the second segment utilized/ 2 = 0.07. To 
accomplish this we used equation (10), using the variable values found within the 
paper. Species i = 1, ... m was designated as H + ions only, while /' — m + 1, ... / 
included all other species in solution. 

Calculations conducted on bare, insulating piezoelectric materials used the same 
methodology as those done for metal electrodes-covered insulating piezoelectric 
materials. The only difference is that the valence band and conduction band energies 
of the piezoelectric material were used in place of the metal's Fermi energy (equation 
16 and Table 1). In the case of the semiconducting piezoelectric, the same procedure 
was followed the case of the insulating piezoelectric [equation (18)] except that V Semi 
[equation (17)) was used in place of V P max . 

The experimental data was acquired by using a PMN-PT single crystal slab as the 
piezoelectric component for piezocatalyzed water splitting (from our previously 
publication 1 ). The piezoelectric cantilever architecture was constructed and placed 
within a sealed chamber with access ports for piezoelectric actuation and monitoring, 
environmental purging, and atmospheric sampling. Straining of the piezoelectric 
cantilever was achieved by external electronic actuators. The H 2 concentration was 
determined by an AMETEK Analyzer ta3000F H 2 gas analyzer. The piezoelectric 
potential was measured by a digital oscilloscope and a potentiostat, respectively. See 
reference 1 for more details. 
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